Emergent Run-and- Tumble Behavior in a Simple Model of Chlamydomonas with 

Intrinsic Noise 
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Recent experiments on the green alga Chlamydomonas that swims using synchronized beating 
of a pair of flagella have revealed that it exhibits a run-and-tumble behavior similar to that of 
bacteria such as E. Coli. Using a simple purely hydrodynamic model that incorporates a stroke cycle 
and an intrinsic Gaussian white noise, we show that a stochastic run-and-tumble behavior could 
emerge, due to the nonlinearity of the combined synchronization-rotation-translation dynamics. This 
suggests the intriguing possibility that the alga might exploit nonlinear mechanics — as opposed to 
sophisticated biochemical circuitry as used by bacteria — to control its behavior. 
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Microscopic organisms need to develop swimming 
strategies that can tackle the low Reynolds number con- 
ditions where viscous forces dominate [1]. Moreover, 
their survival crucially depends on developing efficient 
search strategies for the needed chemicals [2[ , as has been 
studied extensively for many bacteria and in particular E. 
Coli [3], which exhibits a run-and-tumble behavior 0, 
The pattern of having relatively long run segments in 
the trajectory intercalated with burstlike tumble events 
during which the orientation of the swimming micro- 
organism is completely randomized is known to have sig- 
nificant advantages as a search strategy jsj. In bacteria, 
the onset of the tumble events is believed to be triggered 
when certain receptors report a change in the concentra- 
tion of specific chemicals, which they measure through 
temporal integration, via signaling pathways that have a 
feedback control on the preferred direction of the rotation 
of the flagellar motors [7j . 

Recent experiments on Chlamydomonas, which is a 
unicellular green alga that swims with two (flexible) flag- 
ella which beat with a breaststroke-like motion [8|, have 
revealed a similar pattern of behavior [9[ . It is observed 
that when the two flagella are synchronized the cell swims 
in a straight line, and that the intervals of synchronous 
swimming are interrupted by periods of asynchronous 
beating that leads to reorientation of the cell. The syn- 
chronous and asynchronous periods of beating are anal- 
ogous to the run-and-tumble motion observed in flagel- 
lated bacteria. This suggests that synchronization of the 
beating flagella — which has been studied using a vari- 
ety of different theoretical and experimental approaches 
over the last decade [9l-[l2| — and its interplay with the 
swimming motion of Chlamydomonas play a central role 
in controlling this stochastic behavior. A key question 
naturally arises: to what extent is a direct biochemical 
switch necessary for the observed run-and-tumble behav- 
ior in Chlamydomonas? In other words, is it possible to 
have an alternative mechanical switch that could lead to 
the same behavior? This is the question we would like to 
address in this Letter. 



We consider a minimal three-sphere model of Chlamy- 
domonas in which one bead models the cell body and 
each flagellum is represented by a bead moving on a tra- 
jectory fixed relative to the cell body, as shown in Fig. 
QJa). The model was inspired by the recent experimental 
measurements of the flow field around Chlamydomonas 
[HI that was found to be represented reasonably well by 
three force monopoles (or Stokeslets) even in the near- 
field regions. Instead of prescribing the motion of the 
beads as in Ref. [14| , we treat the phases that define the 
locations of the flagellar beads along their trajectories as 
dynamical variables as in Ref. and study the com- 
bined dynamics that describes their synchronization (or 
lack thereof) as well as rotational and translational mo- 
tion of the model swimmer. Friedrich and Jiilicher, who 
independently developed the same model in a recent pa- 
per [12[ , showed that the two flagella can synchronize via 
the coupling with the background hydrodynamic flow, 
and that hydrodynamic interaction between the flagella, 
though essential for the net swimming, has a secondary 
role in synchronization. Here, we show that by introduc- 
ing a simple stroke pattern in the cyclic motion of the 
flagellar beads (as opposed to the constant forcing 
considered in Ref. [12(1 ) an d ( a relatively small amount 
of) Gaussian white noise, a run-and-tumble behavior 
emerges from the nonlinear dynamics of the model. We 
observe long run segments in which the model Chlamy- 
domonas swims on straight or slightly curved trajectories 
(depending on the choice of parameters; see below), and 
intermittent sharp tumble events in which its orientation 
is suddenly randomized. For a given noise strength, the 
distribution of the run duration decays exponentially as 
shown in Fig. [lib), in agreement with the observations 
by Polin et al. [9]. The characteristic run time r mn de- 
pends on the noise strength, as shown in Fig. G^c), and 
follows r run oc 1/a, where a is the noise strength. Fig- 
ure [TJd) shows a typical run-and-tumble trajectory for 
completely symmetric driving force (explained below), 
which leads to run segments in which the two flagella 
are synchronized [Fig. (He)] and orientation of the swim- 
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FIG. 1: (color online.) (a) A three-sphere model of Chlamydomonas. The left and right beads represent the flagella and move 
on the circular trajectories shown in the cell frame. The back bead represents the cell body. The faint underlay is a schematic 
of a Chlamydomonas cell, (b) Run duration statistics, in units of cycles, for driving force coefficients ai = a r = 0.7 and noise 
strength a — 10 -4 . The run length distribution decays as exp(— r/r Yun ) with r Tun = 79 cycles in this case, (c) Dependence 
of Trun on noise, for coefficients a\ — a r = 0.7. T run oc 1/cr. The inset shows the fit to a line with slope —1 in log- log scale, 
(d) A run-and-tumble trajectory for ai — a r — 0.7 and a = 10 -4 . (e) 5 versus time corresponding to the trajectory in (d). 
(f) 6 versus time corresponding to the trajectory in (d). The jumps in 6 occur at the same times as the jumps in 5. (g) 
A run-and-tumble trajectory with curved run sections for driving force coefficients a\ — —0.7, a r = 0.7 and noise strength 
a = 0.01. (h) 6 corresponding to the trajectory in (g). (i) 6 corresponding to the trajectory in (g). 



mer is fixed [Fig. DJf)], with large simultaneous jumps 
in the relative phase and the orientation corresponding 
to the intermittent tumble events. A choice of driving 
force with different coefficients (see below) could gener- 
ate curved run segments [Fig. Gig)] via a sequence of 
phase slips [Fig. H{h)], and in conjunction with a simi- 
lar pattern of simultaneous jumps in the phase difference 
[Fig. Hpi)] an d the orientation [Fig. H{i)], in agreement 
with the observations of Polin et al. [9[ . 

The Model. — Consider 3 beads in the x — y plane la- 
beled 'back', 'left', and 'right', which we refer to with 
the subscripts L b\ T, and V, respectively. Each bead 
is of radius a and they are arranged in the configuration 
shown in Fig. [TJa). Let Ro be the origin of the cell frame 
with respect to a lab frame. The cell axes x, y make an 
angle 0(t) with the lab axes X, Y. In the cell frame the 
back bead is fixed; the left and right beads that model 
the flagella move on circular trajectories with radius 6, 
in opposite directions and with phases <pi and <fr r , respec- 



tively. The velocities of the beads are 

R z = Ko + L0y + b(^i-0)ii, (1) 
R r = R - LO y + b(<j) r + 0) t r , (2) 
R 6 = Ko-HOic, (3) 

where the dot denotes differentiation with respect to 
time, and and are unit vectors in the normal and 
tangential directions of the circular trajectory of R^. The 
left and right beads are driven by tangential forces and 
F*, respectively, which define the stroke pattern of the 
cyclic beatings. Normal forces Fp and F™ are exerted 
on the beads in order to constrain them to the circular 
trajectories. The force on the back bead is such that the 
swimmer is force free and torque free: + F r + F5 = 0, 
T/ + T r + T 6 = 0, where F» = F/t^ + Ffni for i = l,r and 
Tj = Hj x Fj for j = 6, Z,r. The forces and velocities 
are related through hydrodynamic interactions between 
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where £ = Girrja is the bead friction coefficient (with n be- 
ing the viscosity of the ambient fluid). In the limit when 
a is considerably smaller than the other length scales, 
the hydrodynamic interaction is described by the Oseen 
tensor 



Gi3 ~ 87r7?|r ij | (I + ^ i ' ii) ' 



(7) 



With Yij = Ti 



The Symmetric Case. — First we consider the case 
where the flagella beat synchronously, and thus <\>\ = 
4> r — (j). We require F{ = F* = F where F = F{(j>) 
can vary periodically with phase and F(<p) > for all 
values of (j). By symmetry, Fp((j>) = F™{4>) — C(0), 
= and R5 = Rbf- First we calculate the constraining 
force C(0), and then solve for </>, which ultimately gives 
us the velocity R^. We find that the cell swims with 
rapid oscillations with a net drift in a direction deter- 
mined by the ratio H/L. We consider the stability of this 
state in the next section. This matches the experimen- 
tal observations made by Racey et al. [16]. The average 
velocity depends on F((j)) and the ratios a/L, H/b and 
L/b. Here, we work with ratios a/L = 1/33, H/b = 1/5 
and L/b = 2. In a real Chlamydomonas cell, a and L 
are of the same order of magnitude, but we choose small 
a/L and a/b in our calculations such that the Oseen ten- 
sor approximation is valid. The velocity scale is set by 
Fo/(6irrja) where Fo = J Q n d(j)F((j))/27r; the kinematics, 
hydrodynamics, and the driving force profile give many 
additional terms to the velocity. The maximum aver- 
age velocity is achieved if the driving force is constant 

Synchronization and Stability. — We introduce a phase 
difference S = <p r — therefore, & = (j) — 5/2 and <p r = 
(j) + 5/2. The evolution of 5 is governed by 




R&) • t r — (Hi — R&) • ti 



+26 cos (5/2) (L cos ^) - H sine 



(8) 



to the leading order. We can solve Eq. (JHj) numerically 
for a choice of stroke pattern (driving force) F* ((/)), i = 
l,r and initial condition 5(<fro) = 5q. Here, we study 
force profiles of the form F*((f)) = Fq(1 + ai cos(/>) where 
-1 < ai < 1. 
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FIG. 2: (color online.) Evolution of 5 for initial condition 
5((j) = 0) = 1.3 with (a) ai = a r = 0.7, (b) a t -0.7, a r 0.7, 
(c) ai = 0.7, a r = 0.8. The insets show a small part of the 
plot in more detail. The bottom axis 0/27T shows the number 
of cycles which increases monotonically with time. 



The behavior depends on the choice of coefficients a\ 
and a r . Figure [2] shows examples of the three main types 
of behavior we observe. When the coefficients are equal 
a\ = a r = a equa i, 5 evolves into the synchronized state 
when 0.5 < a equa i < 0.8 and the initial condition is 
in a region around 5 = 0. Figure [2](a) shows this for 
^equai = 0-7. For other values of a equa i, 5 evolves into 
a periodically oscillating state about n. The synchro- 
nized state is stable for 0.5 < a eqU ai < 0.8 and unstable 
for other values of a eqU ai- We can consider the stabil- 
ity of other driving force profiles with equal coefficients 
by studying how 5 evolves for initial conditions close to 
5 = 0. For example, if we replace the cos^ term in the 
driving force with cos 20, we find that the synchronized 
state is stable for a equa i < —0.3 and unstable otherwise. 

When the coefficients have opposite sign then 5 evolves 
into a periodically oscillating state near 2i\n (where n is 
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an integer). Figure Efb) shows the evolution into a state 
oscillating about 2tt for a\ = — 0.7, a r = 0.7. For other 
values of a\ = — a r , the shape of the oscillation changes 
and the center of the oscillation drifts, but remains close 
to 2irn. When the coefficients have the same sign but 
different magnitudes, then S evolves into a state oscil- 
lating periodically about 7r, as shown in Fig. [2fc) for 
ai 0.7, a r 0.8 jl7|. 

Noise. — The existence of well defined stable dynami- 
cal states suggests that we might obtain sharp stochastic 
transitions between them if we take into account an in- 
trinsic noise in the driving force. We consider the driv- 
ing force as = 1 + (a^ + Ci)cos0i with index 
i = Z,r. d(t) has a Gaussian probability distribution 
with zero mean and correlation function ((i{t)(j(t')) = 
a 2 SijS(t — t'). We find that various types of behavior can 
occur depending on the noise strength a and the choice 
of a/, a r . 

We start by considering the simplest case of a\ = a r = 
0.7, which leads to a quick and robust synchronization 
as seen in Fig. Efa). Figure QJd) and (e) show a cycle 
averaged trajectory of the cell and the evolution of the 
phase difference in time for a = 10 -4 , respectively. The 
phase difference randomly oscillates about 2i\n for a time 
r before slipping away and oscillating about 2i\m\ with 
m ^ n. The distribution of the time between stochastic 
slips r is shown in Fig. QJb) for a = 10 -4 and decays ex- 
ponentially, indicating a Poisson behavior and matching 
the observations of Polin et al. [9] . The characteristic run 
length r run is inversely proportional to the noise strength 
r run oc 1/cr, as seen in Fig. ffl^c). Note that our numerical 
results for larger values of a deviate from the 1 /a behav- 
ior, because for those values we approach the cutoff limit 
in which the run length is comparable to the duration 
of the tumble event. Figure [iff) shows the orientation of 
the cell, demonstrating that the slips in J, shown directly 
above in[]Je), correspond to the changes in orientation. 

We next consider opposite signs on the coefficients, as 
ai = —0.7 and a r — 0.7, which is interesting due to a 
degree of inherent frustration. In the noiseless case, S 
oscillates periodically about 2irn. When noise is added, 
the same behavior is observed as for equal coefficients: 
S oscillates about 2i\n during run phases, then slips and 
oscillates about 2i\m. Figures [TJg), [TJh), andQJi) show 
a trajectory with the corresponding phase difference and 
orientation, respectively, for a = 0.01. Figures []Jg) and 
[TJi) show that the run sections always curve to the left as 
a result of the imbalance, which also gives a bias to slips 
in the positive direction where the right bead performs an 
additional beat. The net swimming velocity is an order of 
magnitude lower for opposite coefficients than for equal 
coefficients. 

We can also obtain run-and-tumble behavior for other 
combinations of the coefficients. For example, for the 
case where we introduce a mismatch in the coefficients, 
say as a\ = 0.7, a r = 0.8 [Fig. EJc); S oscillates about 



(2n + l)7r], the presence of noise introduces phase slips 
and we obtain run-and-tumble behavior where the flag- 
ella beat out of phase during the run segments. 

Discussion. — We obtain run-and-tumble behavior with 
this simple model using only mechanical considerations 
that govern the force-free and torque-free swimming of 
the model Chlamydomonas, dynamical synchronization 
of the two flagella, and intrinsic noise. Remarkably, the 
threshold-like behavior that in bacteria such as E. Coli 
is believed to be controlled by a sophisticated feedback 
mechanism controlled by biochemical signals, emerges 
naturally from nonlinearities of the mechanics of the 
system and could be triggered by uncorrelated white 
noise. The parameters in the driving force do not need 
to be tuned to achieve run-and-tumble behavior, and 
could lead to large variety of behaviors such as straight 
and curved trajectories as well as in-phase and out-of- 
phase synchronization. Another interesting feature of our 
model is that phase slips and persistent rotation in one 
direction could be obtained without the need for an in- 
trinsic frequency difference between the two flagella (i.e 
when the average intrinsic beating frequencies of the two 
flagella are the same but they correspond to opposite 
signs of the coefficients), which could be interpreted as if 
they are going through their strokes with an internal 7r 
phase difference. In other words, the asymmetry between 
the flagella does not necessarily need to be intrinsic and 
could be dynamically enforced. Finally, the same model 
with a slightly different choice of (mismatched) coeffi- 
cients could also describe the case of anti-phase flagellar 
beating during the run segments. This behavior has been 
recently observed in a mutant of Chlamydomonas [181 ]. 

In conclusion, we have presented a simple model for 
Chlamydomonas that can describe all main experimen- 
tally observed features of the flagellar dynamics and 
swimming behavior of the alga. This will hopefully open 
up a wide range of possibilities for quantitative studies 
of their behavior, and help shed light on possibilities to 
exploit mechanical effects and constraints towards bio- 
logical functions. 
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